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Abstract 

ff") The X(3872) has non-charmonium-hke properties, such as decay processes that seem to violate 

<^ isospin, and a mass that hes unexpectedly close to the D^ D^* threshold. An EFT that includes 

(N 

;h both charmonium-like (short distance) and molecule-like (meson bound state) properties is used 

■^ to analyze the A(3872) as it is produced in the decay of V'(4160). This is a route that BESIII 

^rr may be able to measure. We find that the correlation between the angular distribution of the 

I— I outcoming photon (or A(3872)) and the polarization of the '(/'(4160) source may be used to provide 

Qh information on whether short-distance or long-distance effects dominate. 
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I. INTRODUCTION 

The X(3872) was discovered by the Belle collaboration [T] as a narrow resonance from 
the decay B^ — )■ X(3872)i^^, X(3872) — t- J/ip-K^-K^ . Its existence has been confirmed by 
the CDF [2], DO [3], and BaBar [1] collaborations, and now at the LHC [5],|6]. The most 
recent Particle Data Group value for its mass is m(^x) = 3871.68 ± 0.17 [7], but whether 
it is actually above or below the D^ D^* threshold at 3871.81 ± 0.36 MeV is still an open 
question. The Belle collaboration finds an upper limit on the width of the X(3872) to be 
r(x) < 1.2 MeV at a 90 percent confident level 0. 

While uncertain for most of the time since its discovery, the J^*" quantum number as- 
signments for the X(3872) are now known to be l^"*" [ilj. This, along with the closeness of 
the X(3872) to the D^ D^* threshold, makes it possible for the X(3872) to be interpreted as 
a loosely bound state of D^ and D^* mesons. The possibility that mesons could themselves 
form "molecular" bound states of other mesons was discussed in Ref. [10] and for charmed 
mesons in particular in Refs.[TTHIS]- The X(3872) was investigated as a potential molecule 
shortly after its discovery in Refs. P^H2U] . It seems certain that there is at least a component 
of the X(3872) that can be taken as a molecule given that it will likely strongly mix with 
the C=+l combination of the neutral D mesons. Exactly how much of it is molecular, what 
else might describe its wavefunction, and what observables should be studied to unravel it 
are the subject of lively debate in the literature. 

If the X(3872) is indeed a molecule, then it is a very shallow bound state with a very large 
scattering length, possibly in excess of 6 fm. This would make it larger than, for example, 
the deuteron. The benefit of such a shallow bound state is that its properties are dictated 
by this large scattering length. The universal properties of such systems is discussed in 
Ref. [21]. 

In Ref. [22j we explored the behavior of the X(3872) by noting that its production angular 
distribution depends upon the ratio of short-distance to long-distance terms. In particular, 
using the X-EFT developed in Ref. [23] along with heavy hadron chiral perturbation the- 
ory (HHxPT), we determined that the decay of the ?/'(4040) — )■ X(3872)7 depended upon 
diagrams that are dominated by molecular-like (or long-distance) behavior and a single di- 
agram that depends upon a short-distance interaction. Since this is an effective field theory 
(EFT) treatment, whether the short-distance operator mimics a cc character or some other 



short-distance character is not determined. In particular, while the X(3872) may also mix 
with a linear combination of charged D'^*^ mesons as well as neutral ones, we consider those 
"short-distance" (8 MeV above the X(3872) mass) on these scales. 

In this paper we look at production of the X(3872) from the decay '?/'(4160) — )■ X(3872)7. 
BESIII intends to produce and study ^(4160) and in particular use it as a source of X (3872) 
production [24j. Like the '?/'(4040), the '?/'(4160) has quantum numbers J^^ = 1 and is 
likely a traditional charmonium excitation. It is one of the L = D multiplets, 2^Di. Its 
partial fraction to electrons suggest that it may have additional L admixtures, but since 
this is uncertain at the moment we take it to be dominantly a pure state here. Its mass 
and width are estimated by the PDG to be m^ = 4153 ± 3 MeV and T^ = 103 ± 8 MeV 
respectively. 

Below we find the differential cross section -—^ '-^ — — and extract its dependence 

on the angle between the outgoing photon momentum and initial ^'(4160) polarization vector. 
We discuss how this correlation can be used to determine the short-distance versus long- 
distance character of the X(3872). We also provide an estimate for the total decay rate 
r[V^(4160) -> X(3872)7] should the X(3872) be predominantly a bound state of neutral D 
mesons, ^ {D^D*° + D*^D^). 

II. EFT LAGRANGIAN 

To create an effective field theory for QCD, we identify the fields whose behavior we 
want to describe, the energy region of interest, the symmetries we want to impose, and the 
small parameter that will organize the operators in the Lagrangian. Then we write down 
the most general Lagrangian order by order [SH]- In the limit iric — )■ oo and ma^u,s ~^ QCD 
acquires two approximate symmetries: heavy quark spin symmetry and chiral symmetry |2S] ■ 
HHxPT is an effective field theory with both of those symmetries, including a simultaneous 
expansion in both limits. The heavy hadrons are treated as nonrelativistic particles with 
their classical mass term rotated away, leaving a derivative expansion in p/rric, where p is 
the (small) momentum scale in the problem |26j. We will keep the zeroth order terms in 
the chiral expansion, but include the leading p/rric operator. XEFT [23] is an effective field 
theory describing low-energy nonrelativistic D, D, D*, D*, and tt mesons near the D^ + D^* 
mass threshold. It is matched onto HHxPT by integrating out virtual states whose energies 



are widely separated from that threshold. It is similar to the NN-EFT created to treat the 
deuteron as a bound state of nucleons [27], but it is better behaved in that pions can be 
treated perturbatively. XEFT was designed to describe the X(3872) as a bound state of 
D^/D*^ _j_ ^ 1^ mesons. 

For the ^(4160) — > X(3872)7 decay we need HHxPT operators that include the iplAlQO), 
D^*^'^, and Z)*^*)° particles. The D^*^ mesons are collected into a superfield to encode the 
heavy quark symmetry. In the lowest multiplet the quarks cu form a bound state with 
relative orbital angular momentum L = 0. The quark spins combine to form the J = 
D mesons (denoted P) and the J = 1 D* mesons {V^). In general this would include the 
charged and strange-ness containing D mesons as well, but here we only require the neutral 
ones. The superfield is 

i, = l±i(^^„-P-,)i^ , (1) 

where v^ is the heavy quark four-velocity. Because the heavy hadrons are treated as static 
sources, there is no pair production. The Z)(*)° mesons have their own field, H. 

General discussions about combining different spin and orbital angular momentum states 
into one field multiplet can be found in Refs. [2S1ES]- Ref. [2H] provides the multiplet fields 
for the cc states within one L value. Each of the quarks has spin s = 1/2, so the cc state has 
spin S = or S = 1. All the possible J states with the same angular momentum L are then 
given by J = L when S = 0, and J = \L — 1\,L,L + 1 when S = 1. The particle iIj{4160) 
consists of cc quarks which have relative L = 2 angular momentum. So the field multiplet 
in which ?/'(4160) lives is 



^-^{Hrio. + ^{e^'-'^v^^pHl^^ + e'^-^^v^ipH!^;) 






where Ha, K^ are the effective fields of the various members of the multiplet with total spin 
J = A. Since the total spin of ip{4160) is J = 1 we need only the A=l term. 

Using HHxPT power counting we identify the leading order operators that couple D^*' 
mesons to photons, the ■^(4160) to the D^*'^ mesons, and the ip{4:160) to both D^*^ mesons 
and photons. 



£ = ^Tr{H^Ha ■ BQ) + ^Tr{H^a ■ BH) 

(1 — "^ PC — 

+ i-Tr{rH^ai dj H^) + i—Tr{J'UoiEjH) + h.c. , (3) 

where we use the 2-component notation of Ref. [30], with 

H = V -a + P (4) 

the superfield that contains both the vector Vi field of the D*^ and the pseudoscalar field 
P of the D^. Because we are confining ourselves to the neutral D-mesons only, Q=2/3, 
and the isospin subscripts are dropped, aj is the spin Pauli matrix. The second term in 
Eq. ([3]) contains the coupling to the charm quark of Q' = 2/3. For a discussion of these and 
higher order EM couplings among the D-mesons and their excited states see Ref. j3j. In 
the nonrelativistic limit our ■?/'(4160) (the Hi of Eq. d2]))is now called ip^, 



r^ = -J-U^^ + a^^') - ^5'V"^" (5) 

2 V 5 a/Is 



The coefficients /3, g, and c will be discussed in the "Estimates of Parameters" section. 

III. THE DECAY ^'(4160) -^ X{3872)-f 

The Feynman diagrams contributing to the decay are given in Fig. 1. Using the rules 
obtained from Eq. pi) we find each contributes the following amplitude: ^ 

/ N gP+e 1 



3vl5-c/'y — A 
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(^)^-/|^^(A^-e.)(e.^-^^xe;) (7) 



(d) = -o\/ oCcE^e'D* ■ e^ ^ S* , (9) 



^ These were also calculated by T. Melien (unpublished). 
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FIG. 1: Feynman diagrams contributing to the decay ^^(4160) — )• D^D*^ . The thick sohd hne is 
the V'(4160) particle, the thin sohd hne is a D^ particle, the double line is a D^* particle and the 
wavy line is a photon. 



where I3± = (3 ± ^; the polarization vectors of the photon, D"*, and ■?/'(4160) are ey, eo* 
and e*^, respectively; and k is the outgoing photon momentum. The decay rate depending 
on the polarization of the initial ?/'(4160) is found by summing over the final photon and 
D^* particle polarizations: 



r(e-V) ~ ^(^ + Cf\k ■ e-^p + ^-{B - Cf\k x e^p 
where k is the unit vector in the direction of the photon's 3-momentum, and 



A = gP^e^E' 



7^7 - 3A 



B 



3^ ^ A2 - i?2 ' 

ge -P+E^ + 5^.E^{A-E^ 
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Averaging this over iIj{4160) polarizations gives the total decay rate 

T^liiA + Cy + iB-Cr) . 



(10) 



(11) 
(12) 
(13) 

(14) 



If the ip{4160) is produced in an electron-positron collider such as BESIII then it is produced 
with a polarization normal to the beam axis in the limit that the electrons can be treated as 
massless helicity eigenstates. So we can use the relationship between the outgoing photon 
(or ultimately X(3872)) momentum k in the '?/'(4160) rest frame with respect to e^p to obtain 
a relationship between k and the beam axis. 
Defining 



and 



P=^(A + Cy (15) 



T=^{B-Cf , (16) 



the angular distribution of the final states is 

d^ . 2n T-2P .^^. 

where 6 is the angle between the photon momentum vector and iIj{4160) polarization. Sub- 
stituting expressions for A, B and C yields 

2 -E-,+5r(i(A-E^) .2 0/ 4 p 7g^-3A x2 

_ U5 A-E-y ^^1 ^l 15-^7 A2-B2 vJ 

P~ (2 -E^+^rp(A-E^) I X2 I 0/ 4 F 7g^-3A ^ ' V '' 

where r^ = /3_//3+ and r/ = ^-. This ratio t] provides a measure of how much of the decay 
behavior and polarization correlation is driven by the long-distance diagrams (Fig. 1(a)- 
1(c)) that depend upon gl3± versus the short- distance contact c-dependent diagram in Fig. 
1(d). The plot of p as a function of A is given in Fig. 2. Varying r^ within reasonable ranges 
does not change the shape of this curve. If the '?/'(4160) — > X (3872)7 decay were driven 
entirely by long-distance physics then r/ = and p ~ —0.8 (for r^ = 1) or p ~ —0.9 (for 
Tjs = 0.66). A measurement of p ~ —1/3 would not be definitive because that is supported 
by either r^ — )■ oo (short-distance dominance) or t] ~ —1.5. But finding p < —0.7 or p > 
would suggest a significant long-distance contribution to the decay. 

XEFT is used to match the ^^(4160) decays in Fig. 1 to V^(4160) -^ X(3872)7, but this 
just provides an overall constant that cancels in the p parameter 

IV. ESTIMATES OF PARAMETERS 

In this section we discuss the parameters 13, g, and c from Eq. ^ as well as the unknown 
matrix element between the X(3872) and the constituents 1/^/2 [D^D*^ + D*°D°). 




FIG. 2: p as a function of parameter r], with r^ = 0.66. 

The coefficient (3 in Eq. (pi) is found from measured electromagnetic decays among the 
D mesons and their excited states. Ref. [30] finds /3~^ = 275 — 375 MeV using conditions 
similar to the ones relevant for this calculation. 

Now we will use some experimental limits and theoretical estimates to provide an order 
of magnitude expectation for the branching fraction expected for ■?/'(4160) — )■ X(3872)7. The 
partial width found in Ref. [22] 

r[X(3872) -^ ij{2S)-f] 



r 



>0.03 



(19) 



tot 



along with the upper limit on the total width given above, provides us with the limit 
r[X(3872) -^ ip{2S)-f] > 0.04 MeV. (But note that Belle does not see this decay [S3].) 
In Ref. [22] we found that 

r[X(3872) ^ ^(2^)7] = 5^|(0|^e.(A)(rP + rP)|X(3872,A))|2 



A 



V2 



.(25) 



^E^m^K25)r , ,. 

367r nix ^ ^ 



(20) 



where we have replaced the A, B, and C coefficients in Ref. [22] with A2, etc., so they won't 



be confused with the A, etc. in Eqs. (11)-(13) above. 



The first term in Eq. (20) is the matrix element that encodes the overlap between the 



X(3872) and the constituents V^ (the vector mesons D*), P (the meson D), etc. This matrix 
element is not known (although it could be estimated using effective range theory were the 
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binding energy well measured) but it appears in all X(3872) production/decay cross sections. 
If we are able to extract it from one measurement, we can then use it in predictions for others. 
The coefficients A2 and B2 depend upon the coupling between the ip{2S) and the D, D*, etc. 
mesons. The coefficient C2 depends upon an unknown short-distance constant. Refs. [Ml 135] 
have estimated the coupling between "ipC^S) and D mesons to be (72 ~ 2 GeV~^/^. (This is 
the same as the (72 of R-^f- I-31], which finds a value as low as 0.55 GeV~^/^.) This g2 is the 
analog of the g coupling in the Lagrangian of Eq. ^ that couples ?/'(4160) to D mesons. If 
we assume that indeed the X(3872) is dominated by a molecular configuration such that we 
can neglect the impact of C2, we find that this along with the experimental limits provide 
a limit on the matrix element squared above to be > 0.09 GeV'^. To estimate the rate of 
ip{4:160) — )■ X(3872)7 still requires that we estimate the g in Eq. ^. This can be attempted 
by comparing the partial widths of the '?/'(4160) to D^*^ + D^*^ estimates in the quark model 
from Ref. [37]. We find that (7 ~ 1 GeV^^^^. Collecting these, and again assuming that the 
X(3872) behavior is dominated by long-distance physics, we can give an order of magnitude 
estimate that the decay rate 

r[V'(4160) -^ X(3872)7] ~ 2 MeV , (21) 

or a branching fraction of about 2%. This would be even larger with the Ref. [3S] value for 
(72- The parameter c might be estimated by saturating with nearby intermediate states, but 
we see that we can learn something about the X(3872) even without specific knowledge of 
c. 

V. SUMMARY 

We have discussed the differential cross section for iIj{A160) — ?■ X(3872)7 by assuming that 
it has a nonzero overlap with a "molecular" bound state 4| (^D^D*^ + D*^D^^ . We argue 
that a measurement of the angular distribution of the photon (or X(3872)) with respect 
to the beam axis can provide information on whether short-distance (charmonium-like) or 
long-distance (molecule-like) behavior dominates in this decay. We have also provided an 
estimate for the branching fraction of the decay route. 
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